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ABSTRACT 

X-ray observations of Seyfert 1 galaxies offer the unique possibility of observing spec- 
ie/)' tral variability on timescales comparable to the dynamical time of the inner accretion 
^ \ flow. They typically show highly variable lightcurves, on a wide range of timescales, 
■ with Power Density Spectra characterized by 'red noise' and a break at low frequen- 
cies. On the other hand, time resolved spectral analysis have established that spectral 
variability on the shortest timescales is important in all these sources, with the spec- 
tra getting softer at high fluxes (in the 2-10 keV band, typically), while the reflection 
component and the iron line often exhibit a complex behaviour. Here we present a 
. model that is able to explain a number of the above mentioned properties in terms of 
• magnetic flares shining above a standard accretion disc and producing the X-ray spec- 
if. \ trum via inverse Compton scattering of soft photons (both intrinsic and reprocessed 
CO ■ thermal emission from the accretion disc and locally produced synchrotron radiation) . 

We show that the fundamental heating event, likely caused by magnetic reconnection, 
must be compact, with typical size comparable to the accretion disc thickness and 
must be triggered at a height at least an order of magnitude larger than its size. The 
fundamental property of our 'thundercloud' model is that the spatial and temporal 
. distribution of flares are not random: the heating of the corona proceed in correlated 

(*S trains of events in an avalanche fashion. The amplitude of the avalanches obeys a 

power-law distribution and determines the size of the active regions where the spec- 
trum is produced. Due to the feedback effect of the X-ray radiation reprocessed in 
the disc, larger active regions produce softer spectra. With our model we simulate 
X-ray lightcurves that reproduce the main observational properties of the Power Den- 
sity Spectra and of the X-ray continuum short-term variability of Seyfert 1 galaxies. 
By comparing them with observations of MGC-6-30-15, we are able to infer that the 
, accretion disc corona in this source must have a large optical depth (tt ii 1-5) and 

5^ ■ small average covering fraction. 

a ~ 
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1 INTRODUCTION 

X-ray observations of accreting black holes, both in Seyfert 
1 galaxies and in galactic sources (GBHC) typically show 
highly variable lightcurves, on a wide range of timescales. 
Time resolved spectral analysis have established that spec- 
tral variability on the shortest timescales is important in 
all these sources. For example, the analysis of average peak 
aligned shots of the GBHC Cygnus X-l (Li, Feng & Chen 
1999) has shown that, during the hard state and the tran- 
sitions between the two states of the source, the spectrum 
evolves from hard to soft to hard again. Such kind of anal- 
ysis, though, is much better suited for Active Galactic Nu- 
clei (AGN). As a matter of fact, the dynamical timescale 
of a Keplerian accretion flow around a black hole, tdyn — 



9 x 1O 3 ot 3/2 (M B h/1O 7 M ) s (where we have introduced the 
dimensionless radius vj = R/Rs = <? R/2GM), is such that 
it is much easier to perform reliable time resolved spec- 
troscopy for the brightest AGN, where a typical flaring event 
may last hours or days. 

There is in fact growing observational evidence for a 
number of Seyfert Is, that during a flare the X-ray spec- 
trum becomes softer as the 2-10 keV flux increases (see e.g. 
Done et al. (2000); Zdziarski & Grand i (2001); Petru cci a t 
al. (2001a); Vaughan & Edelson ( [200 lj ) and Nandra ( |200l| ) 
for a rec ent re view). More recently, Georgantopoulos & Pa- 
padakis ( 2001 ) found similar trends in a number of Seyfert 
2 galaxies observed with RXTE. 

The analysis of the lightcurves in the frequency do- 
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main is usually presented in terms of Power Density Spec- 
tra (PDS). PDS of AGN are in general characterized by 
'red noise' variability (P(f) oc ,f~ 7 , with 1 ^ 7 2) a nd 



a break at low frequencies (Lawrence & Papadakis 1993 



Markowitz & Edelson 2001). Red noise is characteristic of 



man y a strophysical systems (solar flar es (Lu fc Hamilto n 
199l|)7Joft Gamma Repeaters bursts JGogus et al. 199*9^ 



Gamma-ray bursts ( Bcloborodov et al. 200C| ; McBreen et al 
19940~jtc.) as well as more familiar ones (tides, earthquakes 
and even radio broadcasting, see e.g. Press 1978; Montroll & 
Shlesinger 1982), and is intimately related to the turbulent 
and/or chaotic nature of such non-linear systems. Never- 
theless, the information contained in the PDS alone is not 
enough to discriminate between different mechanisms that 
many authors have proposed as the origin of the 1// noise. 

In one of such attempts, Mineshige et al. (1994) (see 
also Mineshige & Negoro 1999) have proposed that black 
holes accretion flows develop into a state of Self-Organized 
Criticality (Bak, Tang & Wiesenfeld 1987) in which mass 
accretion occurs in avalanches, whose size is distributed as 
a power-law: large (and longer) shots contribute to the low- 
frequency part of the PDS, while small and shorter shots 
determine the power-law decline at high frequencies. 

In general, in order to explain the observed values of 
7 in AGN and GBHC, any model requires individual sig- 
nals of very different timescales. In this way Poutanen & 
Fabian ( 1999 ) were able to reproduce the properties of the 
observed Power Spectrum of the GBHC Cygnus X-l using 
a simple stochastic pulse-aval anche model, originally pro- 
posed by Stern and Svensson (1996) to explain 7-ray burst 
temporal characteristics. In such model the individual sig- 
nals correspond to coronal hard X-ray magnetic flares, which 
possess a range of durations and are capable to trigger larger 
avalanches. 

Magnetic field reconnection of flux tubes in a tenuous 
corona is in principle capable of producing power-law dis- 
tribution ofactive region sizes. In pioneering work, Tout & 
Pringle ( 1996 ) have demonstrated both with a simple an- 
alytic model and with Monte Carlo simulation, that, if a 
dynamo-generated magnetic field with a coherence scale of 
the order of the disc thickness emerges from the accretion 
disc, then the inverse cascade process of stochastic magnetic 
reconnection of flux tubes may lead to much larger coherent 
fields and active regions, with a power-law distribution in 
size. 

Recently, Uttley & McHardy (2001) have discovered a 
remarkable linear correlation between RMS variability and 
flux in two X-ray binaries (Cygnus X-l and the millisecond 
pulsar SAX J1808.4-3658). They have also found that the 
lightcurves of three Seyfert galaxies are consistent with such 
relationship and have therefore argued that such behaviour 
is an universal one in compact accreting systems. 

Here we propose a simple model of Seyfert Is high en- 
ergy emission that combines the basic properties the accre- 
tion flows must have to reproduce the short timescales tem- 
poral and spectral characteristics of the observed objects. 

The basic idea of our model is that magnetic reconnec- 
tion in the accretion disc corona does not occur randomly in 
time and space, but rather in avalanches of correlated events 
(Poutanen & Fabian 1999). This not only affects the proper- 
ties of the lightcurve, but, we argue, is responsible for the ob- 
served spectral variability. The spectrum produced in region 



where a large number of heating events take place is differ- 
ent from the spectrum generated by an isolated event, as we 
discuss in the next section. Furthermore, large avalanches of 
magnetic flares may obscure the underlying disc as seen from 
the observer, and complicates the temporal behaviour of the 
secondary spectral features (reflection hump and fluorescent 
lines) . 

To build a full time-dependent model on the basis of 
the above ideas, including a detailed treatment of the ra- 
diative transfer in the accretion disc - corona system is un- 
doubtedly a formidable task. Here we have chosen to sim- 
plify it by adopting an approximate analytic treatment for 
the spectral generation, following the works of Di Matteo, 
Celotti & Fabian (1997; 1999) and of Wardzinski & Zdziarski 
(2000) (as discussed in Appendix A). This, we believe, gives 
a fair representation of the modeled system and enables us 
to highlight the main prediction of the model in the simplest 
possible way. 

The paper is structured as follows: in the next section 
we outline the main features o f ou r model, dealing sepa- 
rately with its spectral (section 2.1) and temporal (section 
2.3) structu re, a nd we elaborate on the thundercloud anal- 
ogy (section 2.2). The results of numerical simulation and a 
discussion of the main physical consequences of the model 
are presented in sections g and ^j, respectively. Finally, in 
section n we draw our conclusions. 



2 MODEL OUTLINE 

2.1 Structured coronae and thermal 
Comptonization 

The best model to date that explains the spectra of Seyfert 
1 galaxies is based on thermal Comptonization of soft pho- 
tons in a hot electron cloud (Shapiro, Lightman & Eardley 
1976; Sunyaev & Titarchuk 1980). Although it was early 
shown that the main features of the high-energy emission 
of black hole powered sources could be explained by the in- 
terplay of the hot cloud (the so-called corona) shining on a 
standard (optically thick, geometrically thin) accretion disc 
(Haardt & Maraschi 1991), the geometry of the hot and 
cold phases is still matter of debate. Such geometry could 
in principle be deduced with sufficient spectral information: 
for example, reflection features produced by a cold, opti- 
cally thick medium illuminated by the hard X-ray power- 
law continuum are sensitive to the detail of the relative 



geometry of the two comp o nents (George fc Fabian 1991 



Matt, Perola fc Piro 1991; Magdziarz fc Zdziarski 1995) 



Simultaneous studies of UV (likely produced by thermal 
emission from the optically thick medium) and X-ray ra- 
diation have already been used in many cases to rule out 
a slab corona in favour of a p atchy one, made of a nurn- 
ber of separate active region s ([Haardt, Maraschi fc Ghis- 



cllini 1997; stern et al. 199E; Di Matteo, Celotti & Fabian 



199E ), while the observed correlation between the photon in- 



dex and the reflection strength (Zdziarski, Lubinski & Smith 



199£ ; Gilfanov, Churazov, fc RcTmvtsev**2*00*o[ ) (which has 



been proved to hold on timescales much longer than the 
dynamical one) has demonstrated the need for a further ge- 
ometrical/dynamical parameter, such as the relativistic bulk 
motion velocity of the coronal material (Beloborodov 1999a; 
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Malzac et al. 2 001) and/or the truncation radius o f the inner 
accretion disc (Zdziarski, Lubinski & Smith 1999). 



As we focus here on short-term variability, we do not 
consider these latter possibilities (namely, we consider a 
static corona and a disc extending down to the innermost 
stable orbit) and assume that the X-ray spectrum is pro- 
duced by thermal Comptonization in spherical active regions 
of size R = rRs, lifted above the disc to a height H = hRs 
(Di Matteo, Celotti & Fabian 1999), and permeated by a 
magnetic field with intensity B. The seed photons for Comp- 
tonization are both synchrotron ones, produced locally by 
the interaction between the hot electrons and the magnetic 
field, and black-body ones coming from the underlying cold 
disc. 

An active region illuminates the cold disc: part of the 
flux is reflected and part is absorbed and reprocessed. The 
geometry of the active region plays a crucial role in deter- 
mining the comptonized spectrum: if L(r) is the luminosity 
of an active region, we have for the reprocessed luminosity 
that returns to the active region and is comptonized there 
(Beloborodov 1999; Malzac, Beloborodov & Poutanen 2001) 



L rep (r) = L(r)(l-a)(l-Mo)/2, 



(1) 



where a is the disc albedo (that will be kept fixed to the 
value of 0.15 throughout the paper), and fxo is a geometric 
factor which regulates the feedback mechanism. As in Di 
Matteo, Celotti & Fabian (1999), we define a circular area 
of radius 7?hot = Hiot-Rs m the accretion disc heated by 
an active region up to a temperature that depends on the 
illuminating flux, and define 

IM> = h/(rl ot + h 2 ) 1/2 (2) 

as the angular size of that circle measured from the center of 
the active region (at height h). As a typical relevant radius 
we take that one at which the temperature has decreased by 
a factor of 2, and we have rhot — 2.3h + r. 

If most of the accretion power is dissipated in the 
corona, and the reprocessed radiation is the dominant source 
of soft photons for Comptonization, the ratio h/r determines 
the amplification factor, A(r), and the active region temper- 
ature (9 = kT c /m c c 2 ): A(r) = L/L Icp = 1 + 46 + 160 2 . 
The resulting photon spectral index, F, can be related to 
A by the approximate expression (Beloborodov 1999a,b) 
F » 2.33(A - I)' 6 , where 5 w 1/6 for GBHC and S » 1/10 
for AGNs. Clearly, for very large active regions (h/r — > 0), 
the amplification factor has a minimum and the spectrum is 
softer. 

In the more general situation the dominant source of 
soft photons can be determined by comparing the local en- 
ergy densities of the synchrotron radiation and of the ther- 
mal emission emerging from the disc, both from intrinsic 
dissipation and reprocessed radiation. As already discussed 
by Di Matteo, Celotti & Fabian (1997, 1999) and Wardzihski 
& Zdziarski (2000), the relevance of synchrotron emission in 
a magnetically dominated active region decreases with in- 
creasing central source mass (so it is less important in AGN) 
and is a strong function of the temperature (the hotter the 
corona the more it is important). Therefore, in the coronal 
flow around an active galactic nucleus, synchrotron radia- 
tion is likely to be relevant only for very compact active 
regions (see Appendix A). 

For each active region size r we calculate the spectrum 



self-consistently: given an active region luminosity L(r), we 
assume that the energy density of the magnetic field is in 
(almost) equipartition with that of the local radiation (Di 
Matteo, Celotti & Fabian 1997): 



8tt 



3eM-L(r)to 
47r(ri? s ) 3 ' 



(3) 



where em is the equipartition factor, of the order of unity, 
and to(R) = R/vdi ss is the dissipation time. The dissipation 
velocity Udiss depends on the uncertain nature of the heat- 
ing process. For magnetic reconnection models it is related 
to the Alfven speed va and to the magnetic Reynolds num- 
ber _R m (which is very large in the low resistivity coronal 
plasma) and ra nges from va/VR™, in Sw eet-Parker recon- 
nection models (Sweet 1958; Parker 1979|), to y x/hiRm, in 
Petschek reconnect ion models (Petschek 1964). In the fol- 
lowing (see section 2.3), on the basis of time variability ar- 
gument, we will show that i>di ss cannot be larger than a few 
percent of the speed of light, and discuss further how to 
constrain the dissipation time. 

We calculate the self-absorbed synchrotron luminosity 
using the expression given in Wardzihski & Zdziarski (2000); 
the intrinsic disc luminosity, instead, is simply given by 
Lmt = rh(l — fnjL-Edd, where rh is the accretion rate in 
units of the Eddington one, /h is the fraction of the accre- 
tion power dissipated in the corona (Svensson & Zdziarski 
1994) and L Ed d = 1.3 x W 3S M/M Q is the Eddington lumi- 
nosity. Finally, the reprocessed luminosity, L Tcp (r), is cal- 
culated from Eq. (|l|). The soft photon input flux from the 
accretion disc is then given by the sum of the two thermal 
components at different temperatures, the intrinsic and the 
reprocessed one (see Appendix A for details). 

In the following calculations of our fiducial model we 
will fix the overall normalization constants Mbh = 10 7 Mq, 
rh — 0.1 and /h = 0.7. Furthermore, we neglect radial de- 
pendence of the energy output and consider only the coronal 
flow above the inner, more luminous part of the accretion 
disc, so that -Rdisc = 20i?s- 

To calculate the comptonized spectrum we use the an- 
alytic expression given in Wardzihski & Zdziarski (2000) for 
a homogeneous and isotropic source characterized by optical 
depth tt and temperature 0. To derive the emerging spec- 
trum we fix the coronal optical depth, t t, th e active region 
size r and its luminosity L(r) (see section 2.2), and calculate 



the temperature by imposing that the total luminosity of the 
active region (synchrotron plus inverse Compton emission) 
equals L(r) (as discussed in Appendix A). 

In Fig. 1 we plot the coronal temperature, 0, and the 
slope of the comptonized continuum, a, emerging from each 
active region as functions of an active region size, for differ- 
ent values of the optical depth. As anticipated, larger and 
more luminous regions are cooler, and henceforth softer, due 
to the larger reprocessed soft radiation they intercept. Al- 
though calculated for a spherical active region, for the sake 
of consistency we have checked that, in the limits h/r — > oo 
and /h — > 1, our results agree with the prediction of Haardt 
& Maraschi (1993) for a slab corona in energy equilibrium 
above a cold accretion disc. 
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Figure 1. The equilibrium temperature (© = kT/m e c 2 , solid 
lines) and spectral index a = T — 1 (dashed lines) of an active 
region as a function of its size, for a fixed height H = lijg and 
different opt ical depths, calculated for C = 10~ 3 and D = 1 
(see section |2.2| ). fn order of decreasing line-thickness, the dif- 
ferent curves correspond to tt = 0.5,1,2. Due to the increased 
reprocessed radiation that is intercepted in the larger regions, the 
temperature drops with increasing r, and the emerging spectrum 
is softer. The shoulder at low values of R for the tt = 0.5 case 
is due to the increased importance of synchrotron emission as a 
source of soft photons (see also Figure Al). 



2.2 Active regions distribution and luminosity: 
the thundercloud analogy 

In the hot, comptonizing corona the energy input is likely to 
be provided by a strong, highly intermittent magnetic field 
(see Merloni & Fabian 2001) amplified in the underlying 
turbulent disc by a dynamo process generated by the Mag- 
neto Rotational Instability (Balbus & Hawley 1998; Tout 
& Pringle 1992). Dissipation is then triggered by magnetic 
reconnection of flux tubes on the smallest possible scale 
(given approximately by the disk scaleheight) . If these small 
reconnection sites were distributed randomly in time and 
space, all with a similar aspect ratio, h/r, we would observe 
shot-noise-like variability in the lightcurve and, more im- 
portantly, no spectral variability. If, on the other hand, the 
system is in some kind of critical state, such that each small 
flare can trigger an avalanche in its immediate neighborhood 
(i.e. the flares are spatially and temporally correlated), the 
situation will be qualitatively different. In fact, any such 
avalanche is effectively a large active region: more luminous 
(because it contains more reconnection sites), more extended 
and with a softer spectrum (because h/r is smaller, see pre- 
vious section). 

As a consequence of such avalanche generation process 
at any given time there would be a spread of active regions 
above the disc, with typical size ranging from the very small 
micro-flares (R ~ disc scaleheight), to the larger avalanches 
(R ~ a few Schwarzschild radii) . 

Hence, we assume that the luminosity of an active re- 
gion is determined by the size of the avalanche, and therefore 
scale as the size of the region itself 

L(r) = cim/H£Edd»- D . (4) 
Here D may be related to the spatial distribution of the 



correlated flares (and be regarded as a fractal dimension of 
the reconnection sites) and/or to the radial dependence of 
the energy generation in the accretion disc. 

On the other hand, because of the stochastic nature of 
the avalanche generation process, we will assume a power- 
law shape for the time-independent active region size distri- 
bution n(r) (as also suggested by the shape of the observed 
power density spectra), such that the number of active re- 
gions of size r lying in the interval (r, r + dr) is, at any time, 
given by 

n(r)dr — C2r~ p dr. (5) 

The two constants ci and C2 can be fixed imposing the over- 
all normalization for the total corona average covering frac- 
tion C and for the average hard luminosity: 

C = n^n(r)^±dr (6) 

J r . ^disc 

f r max 

ihard = m/H-t/Edd = / n(r)L(r)dr, (7) 

J r min 

where A(r) is the area of the disc covered by an active region 
of size r and Adi sc is the total area of the disc. 

To summarize, we model the accretion disc corona as 
highly inhomogeneous stochastic system, whose basic build- 
ing blocks, the active regions, can be viewed as 'magnetic 
thunderclouds', charged by the differential rotation of the 
underlying disc and/or the turbulent motions in the accre- 
tion flow. The sizes of the thunderclouds are distributed as a 
power-law. The fast energy release, triggered by magnetic re- 
connection on the smallest scales, heats progressively larger 
active regions. Each active region (thundercloud) produces 
the observed rapid flares (X-ray lightning strokes) by in- 
verse Compton scattering soft photons coming mainly from 
the underlying optically thick accretion discQ. Furthermore, 
if the coronal optical depth t t is high enough, the active re- 
gions may obscure the X-ray spectral features produced in 
the cold disc (in particular the fluorescent lines and the re- 
flection hump, as discussed in Matt, Fabian & Reynolds 1997 
and Petrucci et al. 2001b, respectively) for an observer sit- 
uated above them (and act therefore as 'Compton clouds'). 

2.3 Time variability and Power Density Spectra 

In general, if the heating process is related to the dissipation 
of the magnetic field in the corona, the dissipation timescale 
is given by t (R) ~ 10 2 rM 7 b s, where M 7 = M BH /1O 7 M 
and b = c/vdis (Haardt, Maraschi & Ghisellini 1994). Here 
we assume that the spectrum is produced instantaneously in 

* It is interesting here to remark that the analogy between high 
energy astrophysical processes and terrestrial lightning is not new. 
McBreen et al. (1994) have pointed out the similarity between the 
log-normal properties of gamma-ray bursts lightcurves and those 
of terrestrial lightning. In the latter case, the durations, peak cur- 
rents, intervals between the strokes in the flashes, and the flash 
charges are log-normally distributed (Berger et al. 1975). When 
the dispersion of a log-normal distribution is large, the distri- 
bution is mimicked by a i/i distribution over a wide range of x, 
while it can be shown that amplification processes generally relate 
log-normally distributed variables to power-law tails (Montroll & 
Shlesinger 1982). 
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every active region, i.e. the timescale over which an active 
region is heated is longer than the region ligh t- crossing time 

for a discussion 



30 (see section 4.2 



and b a few. We fix b 
of this point). 

Power density spectra of both AGNs and GBHCs show 
a break that separates low-frequency white noise (P(f) ~ 
constant) from the red noise part of the PDS. Within our 
model, such a break can naturally be associated with the 
largest possible active region, and thus fix r max - On the other 
hand, in any stochastic process in which the characteristic 
durations and amplitudes are distributed as power-laws, we 
have 

N(t )I 2 (t ) oc t Q for t ,min < to < to . max (8) 

where N(to) is the probability and I 2 (to) the mean square 
amplitude of a flare with a duration to and the PDS is ap- 
proximately given by 



P{f)<*n 7 = (3 + /3), 



(9) 



in the range (27rto,max) < / < (2-7rfo,min) (Lochner, Swank & 
Szymkowiak 1991). Then, according to our model, for which 
N(to) oc tQ P and I(to) oc t Q , the observed slope of the PDS 
determines the relation between p and D: 



p = 2D + 3-7. 



(10) 



In the following section we present a numerical simulation 
of the thundercloud model for AGN spectral/temporal vari- 
ability, keeping fixed 7 = 1.5 (as usually dictated by the 
observations of AGN), and we discuss the main observable 
consequences of the thundercloud model. 



3 THE FULL PICTURE: RESULTS FROM 
SIMULATED LIGHTCURVES 

In this section we describe the results we obtained from 
simulated lightcurves with variable X-ray spectra. To take 
into account time dependence we introduce a new distri- 
bution n(r,t), defined in such a way that the number of 
active regions of size in the interval (r, r + dr) generated 
in the time bin (t, t + dt) is given by n(r, t) dr dt. In a sta- 
tionary state n(r, t) does not depend on t, and we assume 
n(r) = n(r,t)to(r), where to(r) is the total duration of an 
avalanche. 

A simulated lightcurve (see Fig. Q for an example) is 
then produced requiring that in any time bin the actual 
number J\T(r, t) of avalanches generated with size in the in- 
terval (r, r + dr) be picked from a Poissonian distribution 
with mean n(r, t)dr and then summing over all the active 
regions generated at any time. For the sake of simplicity we 
adopted a triangular shape for the individual flares profile. 
Because we enforce relation ( |lC)[ ) for the indexes of the distri- 
bution, the final PDS has the correct slope in the red-noise 
part (see Figure ^). 

The emerging spectrum is calculated superimposing the 
contribution from all the flaring regions active at that time. 
The photon index V is calculated for every time bin with a 
least square fit to the spectrum with a power-law model in 
the 3-10 keV band. 

We have fixed r m i n = 0.02 and r max = 4. The average 
covering factor C, in turn, fixes the total number of active 
regions and the degree of variability (RMS) of the PDS. 




2x10* 3x10* 
Time (s) 



4x10* 



5x10* 



Figure 2. A simulated lightcurve (bottom panel) with its cor- 
responding evolution of the covering fraction (top panel), for a 
Mbh = 10 7 Mq central black hole accreting at one tenth of the 
Eddington rate. 
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Figure 3. Power density spectrum of the simulated lightcurve 
shown in Fig. ^| The indexes of the active region distributions 
are D = 1 and p = 3.5. To reduce the scatter the PDS has been 
obtained averaging the power density spectra of 16 lightcurves of 
the total duration of 1.2 X 10 s seconds each. 



To explore the parameter space, we have simulated 
AGN light curves of total duration equal to 2.4 x 10 6 s, 
with a temporal resolution of 60 s. We varied the coronal 
optical depth (in the range 0.5 < tt < 2), the average cov- 
ering fraction (5 x 10 -4 < C < 2 x 10~ 3 ) and the luminosity 
scaling index (0.6 < D < 1.4). 
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3.1 The spectral index-luminosity correlation 

In Figure ^, for D = 1 and C = 10~ 3 , we plot the photon in- 
dex r of the total emerging spectrum, measured in the 3-10 
keV band, as a function of the luminosity in the same band, 
for three different values of the optical depth. As discussed 
in section |^, the higher the 3-10 keV luminosity, the softer 
the spectrum in that energy band, as observed. Furthermore, 
when the largest flares occur, and consequently the instan- 
taneous covering fraction has a maximum (see Fig. |), the 
r — L relation tends to saturate to the softest value corre- 
sponding to an almost slab-like geometry. This is an intrin- 
sic property of our model, and is due to the fact that larger 
active regions are cooler because the reprocessed soft lumi- 
nosity that reenters the blobs (and is comptonized there) is 
larger. 

Overall, we find a spectral variability range Ar ~ 
0.3 — 0.4 for a change in luminosity of about a factor of 
four, similar to what is observed in many Seyfert Is. It is 
worth remarking that a similar correlation was indeed pre- 
dicted by Haardt, Maraschi & Ghisellini (1997) for the case 
of a pair dominated corona, but in that case the same AF 
was achieved for a luminosity variation of at least two orders 
of magnitude, which is clearly inconsistent with the obser- 
vations. To verify that, we show in Figure ^| the observed 
r — L correlation for the Seyfert 1 galaxy MCG-6-30-15 
(from Vaughan & Edelson 2001) observed with RXTE in 

19970. 

For both our simulated lightcurves and the MCG-6 data 
we have fitted the F-luminosity relation with the function 

V = F -KL S . (11) 

In Table 1 the values of the three fitting parameters To, K 
and 8 are listed for our set of simulated lightcurves. The 
asymptotic value of the photon index, To, depends mainly 
on the optical depth, and is larger for larger values of Tr. The 
exponent 8, instead, which determines the amount of spec- 
tral variation for a given increase in luminosity, is mainly 
dependent on the index D. 

For the observed data plotted in Figure we obtain 
To = 2.30±g:||, K = 6.8lf 5 8 and 8 = - 1.311 ° 7 7 9 9 . Com- 
pared to our model simulations, such a result is a strong 
indication of a relatively high value for the coronal optical 
depth (tt i; 1.5) in this source, at least for the observa- 
tion considered here. This, in turn, agrees well with the re- 
sults presented in Guainazzi et al. (1999) of the analysis of 
a BeppoSAX observation of the same source: the measured 
spectral cut-off energy (E c ~ 130 ± 40 keV for their best- 
fitting model) , in the framework of thermal Comptonization 
models implies a coronal temperature^ O ~ 0.1, that we ob- 

t A thorough comparison of the thundercloud model with ASCA 
data for MCG— 6-30-15 will be presented elsewhere (Shih et al. 
2001). 

t We note that, usually, an e-folding power-law is not a good 
approximation for the actual shape of the observed high energy 
cut-off in many sources, in particular those with the better quality 
data, see e.g. Frontera et al. (2001). The actual coronal temper- 
ature is somewhat model and geometry dependent, and its real 
value for MCG-6-30-15 may differ from the one we adopt here by 
a factor of the order of unity, which does not change our main 
conclusion regarding the source optical depth. 
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Figure 4. The observed photon index in the 3-10 keV band plot- 
ted versus the luminosity in the same band (in units of 0.1-Le<m) 
as obtained from our fiducial simulation for D = 1 and C = 10 -3 . 
The points for all the time intervals have been grouped in lumi- 
nosity bins and the average gamma has been calculated. The er- 
ror bars correspond to the weighted statistical uncertainty. Also 
shown are the best fitting functions (as in Table 1) for each of 
the four groups. Triangles corresponds to tt = 0.5, squares to 
tt = 1, circles to tt = 1.5 and stars to tt = 2. 

tain for tt ii 1.5 (see Fig. [l]). Although it is more difficult 
to put tight constraints on the other two model parameters 
(D and C), nonetheless we may conclude that the observed 
r — L correlation is consistent with D ~ 1 and C ~ 10 -3 . 
Such a fairly small value for the average covering fraction is 
indeed required by the observed data: the smaller the cov- 
ering fraction, the larger the observer variability, and the 
greater the chance of a large flare to occur. Too a large value 
of C would therefore drastically reduce the observed spectral 
variability. For our fiducial case, D = 1 and C = 10~ 3 , we 
have ci = 0.16 and C2 = 0.046. The total number of active 
regions at any time is on average^] Af= j n(r)dr ~ 160. 

3.2 The variability-luminosity correlation 

As already mentioned in the introduction, recent analysis of 
time variability properties of X-ray binaries and AGN (Utt- 
ley & McHardy 2001) have highlighted an intrinsic property 
of broadband noise variability in accerting compact objects, 
namely the strong linear correlation between RMS variabil- 
ity and flux. 

We have analysed our simulated lightcurves to look for 
a similar correlation. With the parameters D and C fixed to 
the vaues of 1 and 0.001, respectively, we have simulated a 
lightcurve of the total duration of 4.3 x 10 6 s, subdivided it 
in 560 segments of 7680 s each in the first case, and in 2240 
segments of 1920 s each in the second case, and calculated 
the mean luminosity and the variance a for each segment in 
the standard way. The results, binnned in luminosity inter- 
vals, are shown in Fig. |(| A clear correlation is evident in 

S The total number of individual micro-flares that are the funda- 
mental constituents of an avalanche is of course larger, depending 
on the actual number of events that are contained in any larger 
active region. 



© 0000 RAS, MNRAS 000, 000-000 



Spectral and temporal variability of Seyfert Is 7 



Table 1. Results from simulations and best fit parameters for 
the T-luminosity relation. We fitted a function of the form T = 

r - kl s . 
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Figure 5. The observed photon index in the 2-10 keV band plot- 
ted versus the RXTE count rate in the same band for the Seyfert 
1 galaxy MCG-6-30-15, from the analysis of Vaughan & Edel- 
son (2001) (their Figure 3). The data point have been grouped 
in flux bins and the errors correspond to the weighted statisti- 
cal uncertainty on the value of the spectral index. Also shown is 
the best fitting function of the form T = To — KF S . We obtain 
To = 2.30+.2? 3 , K = 6.8±l° 5 a and 8 = -1.31±g;£. 



both cases, consistent with a linear one. We fitted a function 
of the form a = fc(Z/3_io — Lq), with k and Lq constants, 
and found k = 1.49 ± 0.05 and L = 0.091 ± 0.006 for the 
long interval case (dashed line in Fig. []), and k = 1.11±0.04 
and Lo = 0.063 ±0.005 for the short interval case (solid line 
in Fig. []). In the framework of the thundercloud model such 
behaviour can be explained by the non-random nature of 
the flare distribution. In fact, the large number of small, 
hard, low luminosity active regions consitutes a roughly 
constant background, with very low RMS variability, while 
much higher variability is associated with the rare large lu- 
minous events. The different slopes of the correlation in the 
two cases simply reflect the fact that the length of the seg- 
ments in which we subdivided the lightcurve is shorter than 
the longest flare timescale, and we are sampling the PDS in 
its red-noise part. 



4 DISCUSSION 

4.1 Coronal scaleheight 

In the above calculation we have kept fixed the value of the 
height of the flares above the accretion disc. In fact, as long 
as the main source of soft photons for Comptonization is the 
reprocessed radiation from the disc itself (as in the case of 
our fiducial model for which c r > a 3> c s , see Appendix A), 
what matters for the spectral properties of the source is the 
ratio h/r of an active region height over its size. Then, in 
order for the model to accommodate spectral variations as 
large as those observed in many Seyfert 1 galaxies, we need 

' mm ' max ■ 

If the small-scale flux loops are generated by a self- 
sustaining magnetic dynamo operating inside the disc, it is 
natural to identify the size of the smallest active regions 
emerging from the disc, i.e. the micro-flares that trigger the 
larger avalanches, with the accretion disc thickness (Tout 
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Figure 6. Luminosity dependence of the variance calculted for 
segments of a simulated lightcurve (D = 1, C = 10 — 3 and rx = 1) 
7.7x 10 s (filled circles, dashed best fitting line) and 1.9x fO 3 (open 
squares, solid best fitting line) seconds long, with resolution of 60 
seconds. The fitting functions are of the form a = fc(L3_io — Lq), 
with k and Lq constants. 



Prin gfe 1996; Galeev, Rosner fc Vaiana 197£ ; Nayakshin 



& Kaz;,nas 2001). Thus, the observed spectral variability 
implies that the average elevation of a reconnection site is 
larger (and perhaps much larger) than the disc scaleheight, 
as firstly proposed by Di Matteo (1998). 

It is interesting that full MHD simulations of magneto 
rotational instability in weakly magnetized accretion discs 
indeed show the formation of a strongly magnetized corona, 
considerably more extended in the vertical direction than 
the disc itself (Miller & Stone 2000). Moreover, in the frame- 
work of a model similar to the one we present here, Di Mat- 
teo, Celotti & Fabian (1999) studied observations of the 
GBHC GX 339-4 in different spectral states to place lim- 
its on the model parameters, and in particular on the flare 
height. They conclude that the X-ray emission from galac- 
tic black holes in their hard state is more likely a result of 
inverse Compton scattering of soft photons in hot magnetic 
flares triggered high above the accretion disc. 

The analytic treatment of Svensson & Zdziarski (1994) 
of the accretion disc-corona system allows us to give an 
estimate of the ratio of the coronal to disc scaleheight in 
the inner, radiation pressure dominated, part of the disc. 
From such a calculation, we can conclude that indeed, given 
the observed coronal temperatures, whenever most of the 
accretion power is released in the corona (say, /h <; 0.7), 
the corona is much more extended in the vertical direction 
than the underlying disc. We have 

He '-^ 



£ 16 



9 

1/2 



CT 3 / 2 e 1/2 m _1 JM -1 (i - MY 



(1 - /h) 
0.1 



(12) 



where w is the radial disc co ordin ate in units of 
Schwarzschild radii, J(zu) = 1 — yji/m, and the last in- 
equality is obtained taking the minimum of vj a ^ / ' J(vj). 

In the case of a structured corona the situation ought 
to be more complicated: if the flux tubes emerging from the 
disc carry a magnetic field with strength comparable to the 



equipartition with the internal disc pressure, they will tend 
to expand (reducing the strength of the magnetic field) in 
the less pressurized coronal atmosphere (Parker, 1979, §8.4). 
However, if the flux tubes are twisted, as is likely to be the 
case if they are anchored in the differentially rotating disc 
underneath, their sideways expansion is counterbalanced by 
the generation inside the flux tube itself of an azimuthal field 
component, with a net increase of its total magnetic energy 
content. The flux rope will rise buoyantly (Parker, 1979, 
§9.1, 9.6), keeping a much lower aspect ratio (defined as the 
ratio of the loop length over its diameter). To conclude, given 
the major uncertainties on the magnetic flux tube dynamics 
in the disc-corona interface, we believe that it is not unlikely 
for the height of the magnetic flares in the corona to be 
larger than the disc thickness^ as envisaged in our model, 
and, indeed, needed to explain the large observed spectral 
variability. 



4.2 The magnetic field 

Within the framework of active coronal models for the hard 
X-ray emission of black hole accretion flows, the exact nature 
of the heating process remains the most elusive issue. Here, 
following the argument of Merloni & Fabian (2001), we have 
assumed that the magnetic field is the main repository of the 
energy in the corona and that the active regions are heated 
by the dissipation of such energy, probably via magnetic 
reconnection. 

Although there is a growing consensus on the dynami- 
cal importance of magnetic fields in accretion discs, mainly 
driven by recent progresses in numerical MHD simulation of 
accretion flows and magneto-rotational instability, the de- 
tailed topology of the magnetic field is far beyond the reach 
of our current theoretical understanding. The solar analogy 
has provided theorists with ideas and inspiration, if not with 
quantitative insight, but it is wise to bear in mind that al- 
though we are able to actually image magnetic structures in 
our nearest star, even the heating of the solar corona is still 
not well understood. 

With these caveats in mind, we discuss here some of 
the implications of our stochastic model for the accretion 
disc corona on the strength of the magnetic field and on the 
heating mechanism. 

The ratio of magnetic to thermal energy in an active 
region is given by 



E 

mag 

/E a 



13xl0 3£Mm/H& CirD -i 

T-rO 



(13) 



For the standard set of physical parameters we have 
adopted, as long as the energy in the magnetic field is more 
than a few percent of the equipartition with the emitted 



radiation field, we have E n 



25r D 1 E t h and the active 



It As a final remark, we note that a puzzling characteristic of 
solar coronal loops, so often invoked in analogy with accretion disc 
coronal ones, is their tendency to have a nearly uniform thickness, 
instead of being substantially wider at their tops than at their 
foot-points, as expected from a straightforward application of the 
force-free assumption (Klimchuk, Antiochos & Norton, 2000). 
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regions are indeed magnetically dominated. In fact, the mag- 
netic field intensity can be estimated as 

B ~ 1.6 X IQ^r^-^M; 1 ' 2 (1) G. (14) 

Furthermore, we have va/c ~ 0.03r D ~ 1 6 1/ ' 2 . This translates 
into v A /c ~ 0.1r 2(D - 1)/3 (ln_R m ) 1/3 , in the case of Petschek- 



type reconnection (see section 2.1). It is evident that the slow 
Sweet-Parker reconnection model is not compatible with the 
above scenario, because it would lead to too large a magnetic 
field accumulated in the corona. Our adopted value of the 
parameter b, is instead more appropriate for a relatively fast 
reconnection, as in the Petschek-type model. In principle the 
value of such a parameter could be inferred from temporal 
studies, if we independently knew the physical characteristic 
size of an emitting region, as its value affects primarily the 
timescale of the individual flares. 

To conclude, as already discussed in section ^| we would 
like to remark on the fact that synchrotron radiation is 
hardly important as source of seed photons for Comptoniza- 
tion in AGN, but for the hardest, compact short lived flares. 
However, direct synchrotron light at frequencies near the 
self-absorption cut-off could be observed below the peak of 
the disc quasi-thermal hump in sources where the direct 
emission from the accretion disc is strongly suppressed (i.e. 
for low accretion rates, high fraction of power dissipated in 
the corona and/or high relativistic bulk motion of the corona 
itself). 

4.3 Iron line variability 

The X-ray spectra of accreting black holes, and those of 
Seyfert Is in particular, often show a clear sign of reflected 
radiation in the form of a prominent fluorescent iron Ka line. 
This spectral component is sensitive to the geometry of the 
disc-corona system and can be then used to place more firm 
constraints on it. Often such lines are Doppler and gravita- 
tionally broadened, suggesting that the cold material within 
which they are produced is moving at relativistic speeds in 
the deep potential well of the central black hole (Fabian et 
al. 2000) . In the standard picture in which the reflection fea- 
tures are produced by the accretion disc illuminated by the 
hard X-rays coming from the active regions, the strong ob- 
served continuum variability should drive the line variability 
with a short time-lag. 

This simple picture, though, is challenged by the ob- 
servations, which seem to indicate that, although the iron 
line does indeed vary with the varying flux, its flux is not 
correlated with the changes of the continuum intensity or of 
its slope (Vaughan & Edelson 2001). 

In the framework of the thundercloud model, such com- 
plicated behaviour can be explained by the fact that the 
observed line at any time is the super-position of different 
features produced by active regions illuminating the disc 
with different spectral shape and intensity. The ionization 
properties of the surface of the accretion disc, which in turn 
determine the line characteristics, are strongly dependent on 



the illuminating continuum properties ( 


Nayakshin, Kazanas 


& Kallman 200C ; Done & Nayakshin 20 


01 


Ballantyne, Ross 


& Fabian 2001: 


Ballantyne 2001). Therefore the observed 
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Figure 7. The estimated iron line luminosity, as calculated in 
Eq. as a function of the observed flux in the 3-10 keV band 
(in units of O.lLEdd) f° r the simulation with D = 0.6 and C = 
0.002. The points for all the time intervals have been grouped 
in luminosity bins and the average line flux has been calculated. 
The error bars correspond to the weighted statistical uncertainty. 
Triangles corresponds to tt = 0.5, squares to tt = 1, circles to 
tt = 1.5 and stars to tt = 2. 



Furthermore, as the coronal optical depth becomes 
larger the reflection features are strongly suppressed by 
Compton scattering in the corona itself (Matt et al. 1997; 
Petrucci et al. 2001b). This effect becomes more important 
the larger is the active region. Analogously to what discussed 
in Malzac et al. (2001), we define a geometrical parame- 
ter that determines the fraction of reflected luminosity that 
reenters an active region, fj, s = h/\/ Ar 2 + h 2 . Therefore, we 
can quantify the effect of an obscuring thundercloud by es- 
timating the iron Ka line luminosity produced by any given 
illuminating spectrum of an active region in the following 
way: 



L Ka = (L 7 -2o/2)o[At s + (1 - Ms)(l - fsc(Tr))], 



(15) 



iron line is expected to be complex, and variable on a wide 
range of timescales in a complicated way. 



where L7-20/2 is the luminosity above the iron Ka edge 
emitted towards the disc and P sc is the scattering probabil- 
ity in an active region (see Appendix A). A fraction fi s of the 
reflected X-rays reaches the observed directly, while a frac- 
tion (1 — jJ, s )Psc suffers further Compton scattering in the 
corona and is therefore scattered out of the line, reducing its 
absolute flux. We note here that the above expression may 
be in fact underestimating the amount of cloud obscuration 
because the optical depth tt that enters in the determina- 
tion of P sc is the one calculated from the center of an active 
region and is smaller than that seen by the photons coming 
from the underlying disc. 

In Figure [j] we show the estimated total Fe Ka line lu- 
minosity produced by all the simultaneous active regions for 
a simulation with C = 0.001 and D = 1, plotted as a func- 
tion of the total observed 3-10 keV luminosity. Clearly, in the 
case of high coronal optical depth (tt iJ 1.5), although still 
mildly correlated with the observed continuum at low lumi- 
nosity (small flares and small timescales) , the line variability 
is strongly suppressed at high luminosity (larger flares and 
longer timescales) , to the extent that for a variation of a fac- 
tor of three in the 3-10 keV flux, the line flux may increases 
only by about 50%. 
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In general, even for less dramatic cases, we expect the 
line equivalent width Wkc to decrease with the increasing 
flux. For example, in the simulation presented in Fig. hj 
for r = 1.5, we have Wkc ~ 300eV for L3-10 — 0.1 and 
Wkc ~ 200eV for L3-10 ~ 0.3. Such an high value for 
the unattenuated iron line equivalent width (compared to 
the expected value of ~ 100 eV; George & Fabian 1991) 
can be consistent with a static corona model if the iron is 
overabundant with respect to the solar value, as indeed has 
been shown to be required to exp lain the reflection features 
of MCG-6-30-15 jLee at al. 200q ). 

A full simulation of the temporal behaviour of the iron 
line, which is clearly beyond the scope of this work, will in- 
volve a full treatment of the ionization state of the accretion 
disc atmosphere below different active regions (with differ- 
ent illuminating spectra) and of the obscuring effects of the 
corona itself. Similarly, the effect of Comptonization in the 
corona on the shape and intensity of the reflection hump 
may reveal themselves in time-resolved spectral analysis of 
the hard X-ray emission, though this may p rove a difficult 



task for the existing instruments capabilities ( Petrucci et al 
2OOlb0 



5 CONCLUSION 

We have presented a model for spectral and temporal vari- 
ability of the X-ray emission from accretion disc coronae 
around black holes. We have assumed thermal Comptoniza- 
tion as the main radiation mechanism to produce the ob- 
served X-ray continuum. The cooling radiation comes either 



We have simulated X-ray lightcurves of AGN and stud- 
ied the correlation between the photon index and the X-ray 
luminosity. By comparing our simulations with observation 
of the best studied Seyfert 1 galaxy MGC-6-30-15 we con- 
clude that the model is able to reproduce the observed V- 
luminosity correlation and to put some constraints on the 
model parameter. In particular, to explain the trend ob- 
served in the 1997 RXTE observations analysed by Vaughan 
& Edelson (2001) we need a fairly high coronal optical depth 
(tt iJ 1.5), a small average covering fraction (~ 10 ) and 
a scaling index for the active region luminosity D ~ 1. Fur- 
thermore our model naturally explains the linear correla- 
tion between luminosity and RMS variability (at different 
epochs in the same source), recently discovered by Uttley & 
McHardy (2001) in a number of accretion-powered compact 
objects. 



Our model does not include other important geometri- 
cal or dynamical parameters that may vary on timescales 
much longer than that over which the corona is heated 
(which is close to the dynamical one) and may indeed be 
relevant to determine the average long-term spectral prop- 
erties of these sources. In fact, our model does not repro- 
duce the observed correl ation between the spectral index and 
the reflection features (Zdziarski, Lubinski fc Smith 199!:; 



Lubihski & Zdziarski 2001), and no 'static' corona does, as 



demonstrated by Malzac et al. (2001). If the corona instead 
is in stable relativistic bulk motion (i.e. is the base of a 
jet/outflow), and the outflow velocity changes with time, 
or from source to source, then we expect the average spec- 
tral index to correlate with the amount of reflection (|Bc-| 



underlying accretion disc, or from the reprocessed radiation 
in the cold disc or from synchrotron radiation produced lo- 
cally in the magnetically dominated corona. For the case of 
a typical AGN of 1O 7 M0, accreting at ten percent of the Ed- 
dington ratio and for which most of the power is dissipated 
in the corona (/h <; 0.7), we have shown that the repro- 
cessed radiation is usually the main source of soft photons 
but from the case of more compact and hot active regions, 
where the contribution from synchrotron radiation is impor- 
tant. Nonetheless, the magnetic field in the corona acts as 
the main energy reservoir and dominates the local energy 
budget. 

The basics geometric properties of the coronal flow as- 
sumed in our model and, we believe, required to explain the 
observed spectral and temporal variability of AGN are the 
following: 

• The corona must not be uniform, but structured and 
heated intermittently (Merloni & Fabian 2001); 

• The fundamental heating event, a flare likely caused by 
magnetic reconnection, must be compact, with typical size 
comparable to the accretion disc thickness; 

• The height of the reconnection site must be at least an 
order of magnitude larger than its size; 

• The spatial and temporal distribution of the flares are 
not random, but proceed in correlated trains of events in an 
avalanche fashion; 

• The size of the avalanches determines the size of the 
active regions and their luminosity, and are distributed as a 
power-law. 



from t hje intrinsic dissipation of gravitational energy in the loborodov 19994 |Malzac, Beloborodov fe Poutanen 2"00l| ) 



as observed. Then our model would be applicable as it is only 
to the sources with the softest spectra (and indeed the spec- 
tral parameters of MCG-6-30-15 correspond to negligible 
bulk coronal velocity in the model of Beloborodov 1999a). 



On the other hand, such a correlation has not yet been 
proved to hold on the short timescales we focused our atten- 
tion on. The short-time variability should bear information 
on some fundamental properties of any accreting system in 
which an optically thick and geometrically thin disc coexist 
with a hot comptonizing medium heated by dissipation of a 
structured magnetic field. A dynamic corona, where the re- 
connection sites moves with relativistic velocities and gener- 
ate an outflow whose global properties evolve on timescales 
much longer than the dynamical one, is a straightforward 
generalization of the thundercloud model we have presented 
here, and should be consider in order to better constrain the 
geometry and dynamics of the inner accretion flow around 
black holes. 



Thus, detailed analysis of time-averaged spectra will 
help determining the value of the parameter that is ulti- 
mately responsible for the long-term variability and there- 
fore for the T-R correlation (such as the bulk outflow veloc- 
ity, for example), while spectral and temporal variability on 
short timescales will ultimately probe the physics of coronal 
heating and magnetic reconnection, providing a test for the 
model we have illustrated in this work. 
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APPENDIX A: THE X-RAY SPECTRUM 

We calculate the self-absorbed synchrotron luminosity using 
the expression given in Wardzinski & Zdziarski (2000). The 
intrinsic disc luminosity is simply Li nt = m(l — fu)LEdd, 
where m is the accretion rate in units of the Eddington 
one, /h is the fraction of the accretion power dissipated 
in the corona (Svensson & Zdziarski 1994) and Z/Edd = 
1.3 x 10 38 M/Mq is the Eddington luminosity. Finally, the 
reprocessed luminosity, L Icp (r), is calculated from Eq. (Q). 
For the two corresponding thermal soft photon fluxes from 
the disc Fbb,i and i*bb,r, we assume black-body spectra, with 
temperatures Qi and B r , respectively. They are given by: 

1/4 

©i = ( ) i i (Ai) 



\a B 4itR 2 dis 
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Figure Al. The relative contribution to the cooling luminosity 
from synchrotron radiation (c s , solid lines), intrinsic dissipation 
in the accretion disc (q, dashed lines) and from reprocessed radi- 
ation (c r , dot-dashed lines), as functions of an active region size 
r for D = 1 and a covering fraction of 10~ 3 . The energy densi- 
ties are calculated for a fixed height above the disc H = IRs in 
the case of an AGN with Mbh = 10 7 Mq accreting with rh = 0.1, 
/h = 0.7 and €m = 0.1. The thickness of the lines indicates the op- 
tical depth: in order of decreasing thickness we have tt = 0.5, 1, 2. 
Clearly synchrotron radiation is important only for the more com- 
pact sources with the lowest optical depth, which are hotter and 
produce harder spectra, as shown also in Figure 1. Reprocessed 
radiation energy density, instead, always dominates as a source of 
soft photons for Comptonization for large enough active regions. 
Such behaviour is common for all values of the parameters; as a 
rule, the reprocessed radiation is more important for smaller cov- 
ering fraction and higher values of the fraction of power dissipated 
in the corona /h- 



ume, P sc , depends on the coronal optical depth only, and 
determine the emerging spectral index: 

lnPsc (A5) 



In A 



In spherical geometry (Osterbrock 1974), 



8t| l 



l + exp" 2TT (2rT + l)] 



(A6) 



Finally, we have for the power-law normalization 
2.789i 1 



Np = Anr 



c T F bh 



»,i(©0 (- 



Ci-Fbb, 

2.789, 



, r (e r ) ( 



e 



e 

+ c s F £ 



+ 



(A7) 



where the synchrotron flux F B = 2irm c c 3 Oxt /Ac and x t is 
the turnover energy (in units of m e c 2 ) of the self-absorbed 
synchrotron radiation, that has been calculated as in Di 
Matteo Celotti & Fabian (1997). The factor <j> « (1 + 
40 2 )/(l + 4O0 2 ) approximates the relative normalization 
of the power-law extrapolation and the soft photon input 
peak flux, and is derived form a phenomenological analysis 
(Zdziarski 1986). 

The coronal temperature, and consequently the slope of 
the comptonized continuum emerging from each active re- 
gion, are calculated self-consistently solving for Lc = L(r) = 

Cim/H-t/EddT" 13 . 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^TfjX style file. 



9r 



\m c c 2 J 



L re p (r) 



1/4 



(A2) 



Comparing the energy densities of the three photon 
fields inside each active region, we determine the coefficients 
Cs, Ci and c r , with c s + a + c r = 1, that fix the relative con- 
tribution of synchrotron, intrinsic and reprocessed radiation 
soft photons sources, respectively. This is illustrated in Fig. 



Al for our fiducial model parameters and cm = 0.1. It is 
evident that for r <; 0.1 h the contribution of synchrotron 
radiation to the cooling luminosity is always negligible, in 
particular for the high optical depth (lower temperature) 
cases. 

The luminosity due to the scattered photons can be 
approximated as a sum of a cut-off power-law and Wien 
component (Wardzinski & Zdziarski 2000), 



L c = Ni 



(I)' 



! exp"^ 8 +iV w (| 



exp 



-x/e 



(A3) 



where x — hv/m c c 2 is a dimensionless photon energy, Ap 
and Aw are the normalizations of the power-law and Wien 
tail, respectively, and are related by the ratio (Zdziarski 
1985) 



Aw = r E (a) 
A P r E (2Q + 3) 



p 

J SC 5 



(A4) 



where Te(x) is Euler's Gamma function. 

The scattering probability averaged over the source vol- 
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